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This  paper  describes  a  computational  study  of  the  partial  oxidation  of  ethane  to  ethylene  in  a  short- 
contact-tiine  reactor  (SCTR),  using  a  two-dimensional  computational  fluid  dynamics  model  with  full  heat 
and  mass  transport.  Detailed  heterogeneous  and  homogeneous  chemical  kinetic  mechanisms  are  employed 
to  describe  the  chemistry.  Emphasis  is  placed  on  simulating  recent  experiments  in  which  the  platinum 
catalyst  is  added  to  the  front  face  of  the  reactor  while  it  is  operating  (“on-line”  catalyst  addition).  Our 
simulations  indicate  that  the  fundamental  behavior  of  the  ethane  SCTR  prepared  with  catalyst  added  on¬ 
line  is  the  result  of  coupled  heterogeneous  and  homogeneous  chemical  processes.  It  seems  clear  that  low 
CH4  selectivity  results  from  the  lack  of  heterogeneous  CH4  production  downstream  in  SCTRs  prepared 
with  the  catalyst  added  on-line.  Total  ethane  consumption  and  ethylene  production  rates  are  less  strongly 
affected  because  the  homogeneous  route  for  these  processes  compensates  for  the  loss  of  heterogeneous 
activity,  whereas  the  production  of  methaneis  much  more  effective  as  a  heterogeneous  process.  These 
results  indicate  that  the  improved  performance  observed  as  a  result  of  on-line  catalyst  addition  is  due  to 
a  shift  from  heterogeneous  ethane  decomposition  to  homogeneous  decomposition.  This  limits  the  total 
production  of  methane  while  increasing  the  selectivity  to  ethylene.  In  addition  to  predictions  of  ethane 
conversion  and  ethylene  selectivity,  the  model  also  predicts  the  production  of  all  other  major  products: 
H2O,  Ha,  CH4,  CO,  and  COg. 


Introduction 

Ethylene  production  in  the  United  States  cur¬ 
rently  exceeds  55  billion  pounds  annually,  ranking  it 
among  the  most  important  chemical  commodities. 
Experimental  development  of  so-called  short  con¬ 
tact  time  reactors  (SCTRs)  by  Schmidt  et  al.  has  led 
to  reactor  configurations  capable  of  remarkable  eth¬ 
ane  conversion  and  ethylene  selectivity  [1,4-8].  In 
these  reactors,  mixtures  of  ethane,  oxygen,  and  ni¬ 
trogen  (with  or  without  H2)  flow  through  a  1  cm  long 
reticulated  foam-ceramic  monolith  coated  with  a 
catalytic  metal  such  as  platinum.  Autocatalytic 
operation  is  initiated  by  an  external  heat  source  ap¬ 
plied  temporarily.  The  reactor  yields  a  mixture  of 
ethylene,  CO,  H2,  and  H2O,  with  smaller  amounts 
of  CO2,  CH4,  C2H2,  and  higher  hydrocarbons.  The 
residence  time  in  the  reactor  is  typically  less  than  5 
ms,  and  reactor  exit  temperatures  appear  to  be  in 
the  range  900-950  °C  [1].  Ethane  conversion  and 
ethylene  selectivity  comparable  to  conventional 
steam-cracking  technology  can  be  achieved.  SCTRs 


are  attractive  because  of  their  small  size,  high 
throughput,  and  near  absence  of  coking,  which  com¬ 
monly  plagues  conventional  steam-cracking  reactors. 

A  new  experimental  configuration  identified  by 
Bodke  and  Schmidt  [2]  consists  of  the  addition  of 
catalyst  to  an  operating  reactor,  resulting  in  im¬ 
proved  ethylene  selectivity.  This  so-called  on-line  ad¬ 
dition  of  catalyst  involves  the  dripping  of  very  small 
amounts  of  platinum  salt,  or  other  catalyst  formula¬ 
tions,  directly  onto  the  front  face  of  the  hot  reactor 
monolith.  The  resultant  loading  is  less  than  10%  of 
that  achieved  using  the  conventional  catalyst  prep¬ 
aration  technique,  which  results  in  a  uniformly 
loaded  catalyst.  In  addition  to  improved  perfor¬ 
mance,  this  on-line  process  reduces  the  noble  metal 
requirement  and  permits  rapid  altering  or  regener¬ 
ation  of  the  catalyst. 

In  a  recent  paper,  our  simulations  of  SCTRs  for 
ethane  conversion  to  ethylene  indicated  that  both 
homogeneous  and  heterogeneous  processes  are 
important  for  the  unique  behavior  of  the  SCTR  sys¬ 
tem  [3].  In  this  paper,  we  use  the  same  modeling 
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CATALYTIC  COMBUSTION 


Fig.  1.  Schematics  of  the  experi¬ 
mental  SCTR  showing  the  heat 
shield  and  central  catalytic  section 
(top)  and  the  single-pore  reactor 
model  (bottom).  A  front  heat  shield 
is  not  in  place  during  experiments 
with  on-line  catalyst  addition. 


methodology  to  study  the  on-line  catalyst  reactor 
configuration.  Homogeneous  and  heterogeneous 
chemical-kinetic  models  are  coupled  to  a  full  two- 
dimensional  simulation  of  the  fluid  dynamics  and 
heat  transfer  within  a  single  pore  of  an  SCTR.  We 
use  the  model  to  predict  the  production  of  various 
gas-phase  species  and  to  calculate  ethane  conversion 
and  product  selectivities. 

In  this  paper,  we  focus  on  the  simulation  of  the 
new  experimental  results  in  Ref.  [2]  in  an  attempt 
to  explain  the  fundamental  behavior  of  SCTRs  and 
to  validate  the  chemical  mechanisms  presented  in 
our  previous  work  [3],  We  find  that  the  on-line  ad¬ 
dition  of  catalyst  at  very  low  loading  produces  a  clear 
chemical  effect  related  to  the  catalyst  preparation 
technique.  The  front-face  loading  results  in  a  catalyst 
coverage  that  persists  over  a  distance  shorter  than 
the  1  cm  length  associated  with  uniformly  loaded 
monoliths  prepared  in  advance.  Heterogeneous  eth¬ 
ane  decomposition  is  inhibited,  while  the  rapid  oxi¬ 
dation  of  H2  added  to  the  inlet  feed  is  not.  Hetero¬ 
geneous  oxidation  of  either  ethane  or  hydrogen  is 
responsible  for  heat  release  that  drives  the  endo¬ 
thermic  dehydrogenation  of  ethane  either  in  the  gas 
phase  or  on  the  surface.  Heterogeneous  methane 
production  is  decreased  in  the  on-line  configuration, 
while  the  total  production  of  ethylene  is  not.  These 
findings  are  in  general  agreement  with  the  conclu¬ 
sions  drawn  by  Bodke  and  Schmidt  [2].  In  contrast 
to  the  conclusions  in  Ref  [2],  we  also  found  that  an 
important  contribution  to  the  total  ethylene  produc¬ 
tion  comes  from  homogeneous  dehydrogenation  of 
ethane.  Previous  reports  have  maintained  that  ho¬ 
mogeneous  processes  are  too  slow  to  impact  the 
product  distribution  in  SCTRs  [1],  but  our  results 
indicate  that  both  heterogeneous  and  homogeneous 
processes  are  important  for  ethane  conversion  in 
SCTRs. 

Review  of  Experimental  Investigations 

There  have  been  a  wide  variety  of  experimental 
investigations  of  the  partial  oxidation  of  ethane  to 


ethylene  over  platinum-coated  ceramic-foam  mon¬ 
olithic  supports  in  SCTRs.  Of  particular  interest  is 
the  breakthrough  work  of  Bodke  et  al.  [1]  performed 
prior  to  the  consideration  of  on-line  catalyst  addition 
in  SCTRs.  In  that  paper,  Pt  and  Pt-Sn  alloy  catalysts 
were  used,  and  the  addition  of  H2  to  the  reactant 
feed  mixture  was  reported  for  the  first  time.  An  eth¬ 
ylene  selectivity  of  85%  with  ethane  conversion 
greater  than  70%  was  obtained  using  C2Hg:H2:02 
inlet  mole  ratios  of  2:2:1  over  a  Pt-Sn  catalyst.  The 
authors  concluded  that  the  improved  performance 
was  a  result  of  H2  oxidation  to  H2O  at  the  expense 
of  ethane  oxidation.  This  resulted  in  decreased  se¬ 
lectivity  to  CO  and  CO2.  Heat  released  during  this 
rapid  hydrogen  oxidation  drove  the  conversion  of 
ethane  to  ethylene  with  high  selectivity.  Three  pos¬ 
sible  mechanisms  to  explain  these  observations  were 
proposed  [1]:  (1)  purely  heterogeneous  oxidation 
and  dehydrogenation,  (2)  purely  homogeneous  oxi¬ 
dation  and  dehydrogenation,  and  (3)  heterogeneous 
oxidation  of  hydrogen  followed  by  homogeneous 
ethane  dehydrogenation.  A  purely  homogeneous 
model  [9]  was  employed  by  the  authors  [1]  in  a  non- 
flowing,  transient  analysis  to  address  these  mecha¬ 
nistic  possibilities.  The  most  promising  scenario  was 
heterogeneous  hydrogen  oxidation,  followed  by  ox¬ 
ygen-free  homogeneous  ethane  dehydrogenation. 
However,  the  reported  timescale  (~10  ms)  required 
to  produce  the  measured  conversion  and  selectivity 
by  this  process  is  far  longer  than  that  characteristic 
of  the  SCTR  [1]. 

Because  the  focus  of  this  paper  is  the  simulation 
of  the  experimental  on-line  catalyst  addition  results 
in  Ref.  [2],  we  briefly  summarize  those  results  here. 
Monoliths  prepared  on-line  have  catalyst  applied 
only  over  several  millimeters,  with  10  times  lower 
loading  than  monoliths  prepared  in  advance,  which 
have  catalyst  distributed  over  the  entire  length  of 
the  reactor  (10  mm).  Hydrogen  added  to  a  uni¬ 
formly  loaded  Pt-catalyst  reactor  operating  with  a 
C2Hg:H2:02  feed  ratio  of  2:3:1  results  in  an  ethylene 
selectivity  of  72%,  compared  with  65%  in  the  ab¬ 
sence  of  added  H2.  However,  H2  addition  decreases 
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the  ethane  conversion  from  69%  to  52%  with  a 
C2Hg:H2:02  feed  ratio  of  2:3:1.  When  the  catalyst  is 
added  on-line  to  the  front  portion  of  the  reactor,  the 
ethylene  selectivity  increases  to  nearly  85%  with 
C2Hg:H2:02  =  2:3:1,  while  the  ethane  conversion 
only  decreases  to  55%. 

In  the  simulations  reported  here,  we  show  that  the 
timescales  for  homogeneous  dehydrogenation  in  a 
SCTR  are  indeed  short  enough  to  help  explain  the 
performance  of  these  reactors.  We  present  simula¬ 
tions  that  include  detailed  heterogeneous  and  ho¬ 
mogeneous  chemical-kinetic  processes  coupled  to  a 
hvo-dimensional  flow  field  model  with  full  heat  and 
mass  transport.  The  results  indicate  that  oxidation 
processes  occur  predominantly  on  the  platinum  sur¬ 
face,  with  only  a  minor  oxidative  contribution  from 
the  gas  phase.  In  contrast,  competitive  pathways  for 
ethane  dehydrogenation  and  ethylene  production 
occur  in  the  gas  phase.  In  addition,  we  show  that 
shortening  the  catalytic  section  of  the  reactor  in  the 
simulation  results  in  decreased  heterogeneous  de¬ 
composition  of  ethane  to  form  products  such  as 
CH4.  Shortening  of  the  reactor  in  the  simulation  is 
analogous  to  adding  Pt  catalyst  on-line  in  an  exper¬ 
imental  reactor,  where  the  Pt  salt  solution  dripped 
onto  the  hot  monolith  face  decomposes  immediately 
and  remains  near  the  front  of  the  reactor. 


Methods 

Reactor  Flow  Modeling 

The  basic  experiment  with  which  we  compare  our 
simulations  involves  the  reaction  of  ethane  and  ox¬ 
ygen  with  a  nitrogen  diluent  and  a  hydrogen  additive 
in  an  SCTR,  shown  schematically  in  Fig.  I.  Briefly, 
the  experimental  reactor  consisted  of  an  o-alumina 
(92%  AI2O3,  8%  Si02)  monolith  with  45  pores  per 
inch  (ppi)  coated  with  a  Pt  loading  of  <0.1%  by 
weight  deposited  on-line.  The  total  mass  flow  rate 
into  the  reactor  was  5  slpm,  with  a  reactor  pressure 
of  1.2  atm.  The  nitrogen  mole  fraction  was  30% 
when  there  was  no  hydrogen  in  the  feed,  and  the 
ethane-to-oxygen  mole  fraction  ratio  was  always  2:1. 
Error  bars  are  not  shown  for  the  Bodke  and  Schmidt 
data  [2],  but  it  is  stated  that  the  product  gas-phase 
carbon  and  hydrogen  balances  closed  to  within 
±5%. 

In  this  paper,  we  report  simulations  of  the  ethane 
SCTR  using  full  heat  and  mass  transport  in  a  two- 
dimensional,  axisymmetric  flow  field.  The  fluid  and 
heat  transport  are  modeled  using  FLUENT  [10], 
with  a  similar  methodology  to  that  used  previously 
[3,11,12].  For  this  application,  FLUENT  solves  the 
steady  two-dimensional  Navier-Stokes  equations  for 
laminar  flow  cast  in  a  cylindrical  coordinate  system. 
Calculations  performed  on  a  Silicon  Graphics  Indigo 


II  workstation  required  between  1  and  6  h  to  con¬ 
verge,  depending  on  the  solution  initial  guess. 

The  experimental  reactor  geometry,  along  with  the 
representation  of  a  monolith  pore  used  in  the  cal¬ 
culations,  is  shown  schematically  in  Fig.  1.  The  axi¬ 
symmetric  channel  model  serves  as  a  reasonable 
geometrical  simplification  for  a  single  continuously 
connected  pore  within  the  ceramic-foam  monolith. 
Experiments  involving  alternative  support  geome¬ 
tries  indicate  that  the  physical  form  of  the  support 
does  not  play  a  significant  role  in  the  reactor  perfor¬ 
mance  and  that  straight  channel  monoliths,  although 
more  difficult  to  work  with,  perform  as  well  as  foam 
monoliths  [i3].  In  our  simulations,  the  reactor  mon¬ 
olith  pore  was  I  cm  long  and  could  have  an  active 
catalytic  surface  length  ranging  from  0  to  1  cm.  The 
catalytic  section  was  modeled  as  an  ct-alumina  ce¬ 
ramic  substrate  with  a  monolayer  coverage  of  plati¬ 
num,  resulting  in  a  surface-site  density  of  =  1.64 
X  10^®  sites/ m^  (from  the  density  of  bulk  platinum). 
Of  the  1  cm  length  of  the  catalytic  section  of  the 
reactor,  only  the  first  2  mm  was  considered  to  be 
catalytically  active  in  our  current  simulations.  This 
corresponds  approximately  to  the  size  of  a  drop  of 
catalyst  salt  solution. 

Chemistry  Models 

Using  the  standard  chemical  model  within  FLU¬ 
ENT  [10],  it  is  not  possible  to  calculate  the  fractional 
coverages  of  surface  species.  This  capability  is  nec¬ 
essary  for  incorporating  a  detailed  heterogeneous  re¬ 
action  mechanism.  For  this  reason,  we  coupled  the 
main  executable  FLUENT  code  to  external  FOR¬ 
TRAN  subroutines  [12]  that  model  the  heteroge¬ 
neous  chemistry.  Chemical  reactions  occurring  on 
the  surface  consume  some  gas-phase  species  and 
produce  others.  In  addition,  heat  may  also  be  re¬ 
leased  or  consumed  as  a  result  of  these  reactions. 
The  net  result  of  the  external  subroutines  is  to  con¬ 
struct  source  terms  consistent  with  the  reactions  oc¬ 
curring  on  the  surface  that  are  included  in  the  con¬ 
servation  equations  for  each  of  the  gas-phase  species 
and  for  the  enthalpy.  A  detailed  discussion  concern¬ 
ing  the  formulation  of  source  terms  and  the  deter¬ 
mination  of  surface  coverages  can  be  found  in  Ref 
[3]. 

For  the  gas-phase  chemistry,  we  used  a  method- 
ology  very  similar  to  that  used  for  the  surface  chem¬ 
istry.  A  detailed  kinetic  mechanism  was  employed 
consisting  of  modified  Arrhenius  expressions,  which 
can  be  subject  to  third-body  collision  enhancement 
factors,  low-pressure  treatments,  and  a  Troe  [14]  bi- 
molecular  pressure  fall-off  treatment.  The  gas-phase 
chemistry  capability  within  FLUENT  is  more  ad¬ 
vanced  than  that  for  surface  chemistry;  however,  it 
still  lacks  the  flexibility  to  allow  for  low-pressure  and 
fall-off  treatments.  Therefore,  we  again  relied  on  ex¬ 
ternal  subroutines  to  interpret  the  gas-phase  kinetic 
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Gas  Phase 

C.Ho  C2H4  CHj  CO  CO. 


Surface  Phase  ^ 

C2H4  ±11 

ill  .. 

CiIIj  CHjCH,  CHCH,  CCH3  -«H».  CH,  +  C  CO  <-*  CO, 
±H  ±H  ±H  ±0  ±0 

Fig.  2.  Schematic  representation  of  the  important  sur¬ 
face  reaction  pathways  for  carbon-containing  species.  Ad¬ 
sorption/desorption  processes  are  also  indicated. 

mechanism  and  develop  the  species  and  enthalpy 
source  terms  needed  by  FLUENT. 

The  detailed  homogeneous  and  heterogeneous  ki¬ 
netic  mechanisms  are  identical  to  those  described  in 
our  earlier  paper  [3].  In  summary,  there  are  25  re¬ 
active  gas-phase  species  in  the  homogeneous  mech¬ 
anism,  involved  in  131  reversible  reactions,  with  1 
additional  irreversible  reaction.  This  homogeneous 
mechanism  is  a  subset  of  a  much  larger  mechanism 
developed  by  Marinov  et  al.  [15]  for  the  prediction 
of  combustion  behavior  in  rich  CH4  and  ethane 
flames.  Validation  work  involving  comparison  be¬ 
tween  simulations  using  this  subset  mechanism  and 
additional  experiments  is  the  subject  of  another  ar¬ 
ticle  currently  in  preparation.  The  heterogeneous 
mechanism  consists  of  20  surface  species  and  84  el¬ 
ementary  reactions.  Species  containing  more  than 
two  carbon  atoms  are  not  included  in  either  the  ho¬ 
mogeneous  or  heterogeneous  kinetic  mechanisms, 
and  there  are  no  oxygenated  hydrocarbons  in  the 
heterogeneous  mechanism.  Experiments  measuring 
the  selectivities  to  species  with  more  than  two  car¬ 
bon  atoms  indicated  that  the  total  selectivity  to  these 
species  is  approximately  5%  across  the  entire  range 
of  H2  inlet  feed  ratios  [16].  In  our  simulations,  the 
omission  of  these  species  led  to  an  overprediction  of 
the  ethylene  selectivity  by  approximately  5%.  Oxy¬ 
genated  species  were  not  found  in  experiments  at 
the  conditions  considered  here  [1]  but  were  found 
under  leaner  condition  burning  higher  alkanes  [17]. 

A  schematic  of  the  surface  reaction  pathways  for 
the  major  carbon-containing  species  is  shown  in  Fig. 
2.  Important  features  of  this  new  heterogeneous  ki¬ 
netic  mechanism  include  reversible  ethane  adsorp¬ 
tion  to  form  a  surface  ethane  molecule.  Ethane  un¬ 
dergoes  dehydrogenation  to  form  surface  ethyl 
directly,  and  by  an  oxygen-assisted  route  to  form 
ethyl  and  hydroxyl.  Ethyl  undergoes  further  dehy¬ 
drogenation  to  form  ethylidene.  Ethylene  adsorbs 
reversibly  into  a  re-bound  ethylene  configuration  and 
undergoes  a  reversible  isomerization  to  ethylidene. 
The  primary  route  leading  to  heterogeneous  oxida¬ 
tion  of  carbon  consists  of  ethylidene  dehydrogena¬ 
tion  to  form  ethylidyne.  Ethylidyne  then  undergoes 
carbon-carbon  bond  scission  to  form  methyl  and  sur¬ 
face  carbon.  Surface  carbon  oxidizes  to  form  surface 


CO  that  can  either  desorb  or  undergo  further  oxi¬ 
dation  to  form  surface  CO2.  Surface  methyl  can  re¬ 
combine  with  an  adsorbed  hydrogen  atom,  resulting 
in  CH4  desorption.  Surface  hydrogen  oxidizes  to  hy¬ 
droxyl,  which  can  then  combine  with  another  surface 
hydrogen  to  make  H2O  or  disproportionate  to  form 
H2O  and  O. 

Despite  the  large  number  of  reactions  in  the  het¬ 
erogeneous  mechanism,  there  are  relatively  few  ad¬ 
justable  parameters  required  to  obtain  agreement 
between  the  simulations  and  the  experiments.  Most 
values  for  kinetic  parameters  are  either  taken  di¬ 
rectly  from  previous  work  or  assigned  nominal  values 
which  are  not  adjusted.  A  large  portion  of  the  surface 
carbon  chemistry  is  taken  directly  from  Wolf  et  al. 
[18],  in  which  simulations  were  compared  with  ox¬ 
ygen-free  methane  conversion  to  ethane  experi¬ 
ments.  The  hydrogen  oxidation  subset  was  taken 
without  modification  from  work  we  performed 
which  resulted  in  a  mechanism  explaining  both  ig¬ 
nition  and  steady-state  oxidation  at  high  tempera¬ 
ture.  That  work  is  currently  in  preparation  for  pub¬ 
lication.  The  values  for  the  kinetic  parameters  used 
in  that  mechanism  rely  heavily  on  published  values 
in  the  literature.  Hydrogen  and  O2  adsorption  and 
desoi'ption  are  from  Rinnemo  et  al.  [19],  H2O  de- 
soi'ption  energy  is  from  Fisher  and  Gland  [20],  rad¬ 
ical  species  treatment  is  from  Warnatz  et  al.  [21], 
and  the  surface  oxidation  reaction  energetics  are 
from  Anton  and  Cadogan  [22].  The  carbon  oxidation 
subset  is  very  important  for  the  current  work,  and 
the  pre-exponential  factors  for  surface  carbon  and 
CO  oxidation  and  CO  desorption  are  adjustable  pa¬ 
rameters  in  the  mechanism.  The  energetics  for  the 
carbon  and  CO  oxidation  reactions  are  taken  from 
calculations  by  Shustorovich  and  Sellers  [23]  and  ex¬ 
periments  by  Campbell  et  al.  [24].  The  CO  desorp¬ 
tion  energy  is  taken  from  further  work  by  Campbell 
et  al.  [25].  The  other  adjustable  parameters  are  the 
pre-exponential  factors  for  the  recombination  of 
ethylidene  and  hydrogen  to  form  ethyl  and  the  de¬ 
hydrogenation  of  ethylidene  to  form  ethylidyne. 
These  parameters  influence  the  rate  of  heteroge¬ 
neous  hydrocarbon  decomposition  and  the  produc¬ 
tion  of  CH4  and  CO  through  the  subsequent  ethy¬ 
lidyne  carbon-carbon  bond  scission.  Space 
limitations  prohibit  the  tabular  presentation  of  the 
chemical  mechanisms  developed  and  employed  in 
this  research.  Electronic  versions  of  these  mecha¬ 
nisms  can  be  obtained  from  the  authors  via  e-mail. 

Results  and  Discussion 

Simulation  results  are  presented  in  Figs.  3  and  4 
and  show  that  our  combined  homogeneous/hetero¬ 
geneous  mechanism  can  successfully  reproduce  both 
the  trends  and  absolute  values  of  experimental  data 
reported  in  Ref.  [2].  The  results  in  Fig.  3a  show  se¬ 
lectivity  to  ethylene  increases  for  all  142^02  mixture 
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Eig.  3.  Comparison  between  predicted  and  experimen¬ 
tal  selectivities  and  conversion  as  a  function  of  H2:02  for 
C2Hg:02  =  2.0.  (a)  Ethylene  selectivity,  (b)  methane  se¬ 
lectivity.  and  (c)  ethane  conversion.  Closed  symbols  indi¬ 
cate  experimental  data;  open  connected  symbols  indicate 
simulation  results.  Circles  indicate  catalyst  loaded  uni¬ 
formly;  squares  indicate  on-line  catalyst  addition. 


ratios  when  the  catalyst  is  added  on-line.  This  in¬ 
crease  is  small  when  there  is  no  H2  added  to  the 
feed  and  becomes  larger  as  the  H2  feed  ratio  in¬ 
creases.  Taking  into  account  that  5%  selectivity  to 
>C2  species  is  observed  in  the  experiments,  and  that 
this  selectivity  is  implicitly  included  in  our  calculated 
ethylene  selectivity,  the  quantitative  agreement  is 
quite  good  as  well.  The  predictions  for  the  uniformly 
loaded  catalyst  are  in  excellent  agreement  with  the 
measurements.  The  simulations  of  the  on-line  case 
tend  to  underpredict  the  experimental  value  at  the 
highest  H2  feed  by  approximately  6%  selectivity. 

The  results  in  Fig.  3b  indicate  that  adding  catalyst 
on-line  results  in  a  decrease  in  CH4  selectivity,  which 
is  consistent  with  the  increase  in  ethylene  selectivity 
shown  in  Fig.  3a.  Both  the  simulation  results  and  the 
experiments  [2]  show  that  the  decrease  in  CF14  se¬ 
lectivity  is  more  pronounced  at  the  higher  F12  feed 
ratios.  The  agreement  between  the  predictions  and 
the  experiments  is  again  excellent  when  the  catalyst 
is  loaded  uniformly.  For  the  on-line  case,  the  trend 
is  predicted  quite  well,  with  the  CH4  selectivity  in¬ 
creasing  at  a  much  lower  rate  than  the  uniformly 
loaded  case  as  the  H2  feed  ratio  is  increased.  How¬ 
ever,  we  undciqiredict  the  CH4  selectivity  across  the 
entire  range  of  H2  feed  ratios  by  l%-2%  selectivity. 

The  results  indicate  that  ethane  conversion  in 
SCTRs  is  not  greatly  affected  by  the  on-line  addition 
of  catalyst  (Fig.  3c).  The  experimental  results  [2]  in¬ 
dicate  that  there  is  a  decrease  in  conversion  at  the 
lower  H2  feed  ratios  when  the  catalyst  is  added  on¬ 
line,  but  that  at  the  higher  H2  feeds  there  is  a  slight 
increase.  The  simulations  capture  this  behavior  quite 
well  at  the  lower  H2  feeds  but  do  not  predict  the 
crossover  observed  at  the  higher  feeds.  Neither  the 
physical  mechanism  or  mechanisms  responsible  for 
the  crossover  observed  experimentally  nor  the  rea¬ 
son  or  reasons  for  the  simulations  not  capturing  this 
behavior  are  understood  at  this  time.  Nevertheless, 
the  largest  discrepancy  between  the  predicted  and 
measured  results  is  only  9%  conversion,  or  approx¬ 
imately  15%  relative  error. 

While  the  results  in  Fig.  3  strongly  suggest  that 
the  effect  of  on-line  catalyst  loading  is  to  reduce  the 
opportunity  for  deleterious  (i.e.,  methane-produc¬ 
ing)  heterogeneous  reactions  to  occur  by  minimizing 
the  amount  of  catalyst  in  the  reactor,  a  comparison 
of  the  predicted  gas-phase  and  surface  rates  of  pro¬ 
duction  confirms  this  supposition.  In  Fig.  4,  we  show 
the  calculated  total  molar  rates  of  homogeneous  and 
heterogeneous  consumption  or  production  of  three 
important  gas-phase  species.  Ethane  consumption 
rates  are  displayed  in  Fig.  4a,  which  shows  that  het¬ 
erogeneous  consumption  is  larger  (more  negative) 
than  homogeneous  consumption  for  all  H2  feed  ra¬ 
tios  but  that  competition  between  the  processes  is 
close.  The  ratio  of  heterogeneous  to  homogeneous 
consumption  rates  is  8.5  without  H2  in  the  feed  and 
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Fig.  4.  Calculated  total  molar  reaction  rate  for  homo¬ 
geneous  (closed  symbols)  and  heterogeneous  (open  sym¬ 
bols)  processes.  Uniformly  coated  monolith  results  are 
shown  connected  by  dashed  lines,  and  on-line  results  are 
shown  connected  by  solid  lines,  (a)  Ethane,  (b)  ethylene, 
and  (e)  methane. 


drops  to  ratios  near  unity  at  the  highest  H2  feed  ra¬ 
tios.  The  results  in  Fig.  4a  also  indicate  that  on-line 
catalyst  addition  on  the  front  face  of  the  reactor 
causes  a  larger  fraction  of  the  ethane  to  be  consumed 
homogeneously  and  at  the  expense  of  the  hetero¬ 
geneous  process.  The  shorter  catalytically  active  sec¬ 
tion  associated  with  the  on-line  configuration  inhib¬ 
its  heterogeneous  decomposition  of  ethane,  and  this 
permits  homogeneous  processes  to  compete  more 
effectively  for  the  available  ethane.  Hydrogen-atom 
attack  on  ethane  to  produce  C2H5  and  H2  and  the 
decomposition  of  C2H5  to  form  ethylene  and  H  are 
the  most  important  homogeneous  reactions.  The 
rates  of  these  reactions  are  strongly  influenced  by 
temperature,  and  the  simulations  show  that  the  peak 
temperature  near  the  front  of  the  reactor  is  30  K 
higher  for  the  on-line  case  than  for  the  uniformly 
loaded  case.  This  is  due  in  part  to  the  lack  of  endo¬ 
thermic  hydrocarbon-cracking  surface  reactions  in 
the  downstream  portion  of  the  reactor  in  the  on-line 
configuration. 

The  results  shown  in  Fig.  4b  indicate  that  en¬ 
hanced  homogeneous  decomposition  of  ethane  re¬ 
sults  in  a  greater  homogeneous  production  rate  of 
ethylene.  While  the  heterogeneous  production  rate 
of  ethylene  is  larger  than  the  homogeneous  rate  for 
nearly  all  conditions  studied,  the  on-line  configura¬ 
tion  results  in  a  larger  fraction  of  the  production  oc¬ 
curring  in  the  gas  phase  for  all  H2  feed  ratios.  The 
total  rate  for  ethylene  production  is  slightly  lower 
for  the  on-line  configuration,  yet  the  selectivity  is 
seen  to  increase.  This  is  because  the  production  of 
CH4  is  decreased,  as  we  discuss  next.  In  addition, 
ethane  conversion  has  also  decreased,  and  product 
selectivity  is  based  only  on  the  amount  of  reactant 
converted. 

In  Fig.  4c,  we  compare  calculated  homogeneous 
and  heterogeneous  CH4  production  rates  for  both 
catalyst  preparation  methods.  Of  particular  interest 
is  that  homogeneous  production  of  methane  is  al¬ 
ways  small  and  not  strongly  affected  by  the  prepa¬ 
ration  method.  Conversely,  on-line  preparation 
causes  a  substantial  decrease  in  heterogeneous  CH4 
production  compared  with  uniformly  coated  mono¬ 
liths.  Thus,  the  total  rate  of  CH4  production  is 
sharply  decreased  due  to  on-line  preparation  and  the 
associated  lack  of  catalyst  length  over  which  to  crack 
C2  hydrocarbons  to  methane.  This  occurs  to  the  ex¬ 
tent  that  decreased  selectivity  results  despite  the 
slightly  smaller  amount  of  ethane  converted. 


Conclusions 

Our  simulations  indicate  that  the  fundamental  be¬ 
havior  of  the  ethane  SCTRs  prepared  with  catalyst 
added  on-line  is  the  result  of  coupled  heterogeneous 
and  homogeneous  chemical  processes.  It  seems 
clear  that  low  CH4  selectivity  results  from  the  lack 
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of  heterogeneous  CH4  production  downstream  in 
SCTRs  prepared  with  the  catalyst  added  on-line  to 
the  front  portion  of  the  monolith.  Total  ethane  con¬ 
sumption  and  ethylene  production  rates  are  less 
strongly  affected  because  the  homogeneous  route 
for  these  processes  compensates  for  the  loss  of  het¬ 
erogeneous  activity,  whereas  the  production  of 
methane  is  much  more  effective  as  a  heterogeneous 
process.  These  results  indicate  that  the  improved 
performance  observed  as  a  result  of  on-line  catalyst 
addition  is  due  to  a  shift  from  heterogeneous  ethane 
decomposition  to  homogeneous  decomposition. 
This  limits  the  total  production  of  methane  while 
increasing  the  selectivity  to  ethylene.  If  the  only  ef¬ 
fect  of  on-line  catalyst  addition  was  to  decrease  het¬ 
erogeneous  ethylene  decomposition,  then  an  in¬ 
crease  in  heterogeneous  ethylene  production  would 
result,  which  is  not  what  our  predictions  indicate. 
We  show  a  decrease  in  heterogeneous  ethylene  pro¬ 
duction  compensated  for  by  an  increase  in  homo¬ 
geneous  production.  The  fundamental  chemical  pro¬ 
cesses  that  drive  the  performance  of  ethane  SCTRs 
can  be  summarized  as  heterogeneous  oxidation  oc¬ 
curring  near  the  front  of  the  reactor,  coupled  via  heat 
transport  to  homogeneous  dehydrogenation  of  eth¬ 
ane  to  produce  ethylene  and  H2. 
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COMMENTS 


Don  Hardesty,  Sandia  National  Laboratories,  USA.  In 
your  presentation,  you  mentioned  the  important  role 
played  by  the  heat  transfer  (axial  heat  conduction)  in  the 
process.  Please  elaborate.  Perhaps  a  neglected  point, 
please  elaborate  on  the  importance  of  using  a  fully  elliptic 
solution — ^what  effects/conclusions  would  have  been 
missed  had  you  used  a  parabolic  approach? 

Author’s  Reply.  Accurate  simulation  of  heat  transport  is 
important  because  the  extent  that  homogeneous  reactions 
contribute  to  the  ethane  conversion  depends  strongly  on 
temperature.  A  fully  elliptic  solution  permits  conduction  in 
the  reactor  wall,  so  spatially  separate  regions  where  exo¬ 
thermic  (oxidation)  and  endothermic  (cracking)  reactions 
take  place  can  communicate  thermally.  A  parabolic  solution 
method  would  fail  to  capture  conduction  of  heat  upstream 
of  the  catalytic  zone,  gas  preheating  is  underpredicted,  and 
wall  temperature  at  the  location  of  heat  release  is  overpre¬ 
dicted. 

John  Mantzaras,  Paul  Scherrer  Institute,  Switzerland. 
Could  you  explain  the  reason  for  the  better  agreement  be¬ 
tween  computations  and  experiments  in  the  ethylene  se¬ 


lectivity  and  methane  selectivity  of  the  on-line  catalyst  at 
higher  H2O2  ratios?  Could  you  also  comment  on  the  pres¬ 
ence  of  “hot  spots”  on  the  catalyst  wall  with  increasing  H2 
addition?  As  your  model  includes  heat  conduction  in  the 
solid  substrate  and  the  computations  are  autothermal,  it 
would  be  interesting  to  know  the  level  of  temperature  over¬ 
shoots  in  the  catalyst  with  increasing  H2  addition. 

Author’s  Reply.  As  mentioned  in  the  paper,  we  tend  to 
overpredict  ethylene  selectivity  by  —5%  because  the  gas- 
phase  mechanism  we  employ  does  not  include  C3  or  bigger 
hydrocarbons.  Selectivity  to  these  species  is  —5%  (experi¬ 
mentally).  So,  in  fact,  the  agreement  between  experiments 
and  simulations  is  better  at  lower  H2/O2  ratios  than  at  the 
higher  ratios.  The  reasons  for  this  underprediction  of  eth¬ 
ylene  selectivity  at  higher  H2O2  ratios  in  not  immediately 
clear,  but  we  do  notice  that  selectivity  to  methane  is  over¬ 
predicted  at  these  same  ratios.  Perhaps  our  heterogeneous 
chemical  mechanism  over  prediets  CH3(s)  +  H(s)  — > 
CH4(g)  in  hydrogen-rich  environments.  The  addition  of  H2 
to  the  inlet  feed  induces  a  temperature  overshoot  of  almost 
.30  K. 


